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Abstract 

The classification of the hydrodynamical growth mechanisms for the spherical 

bubbles of the low-temperature phase in cosmological phase transitions is 

completed by showing that the bubbles can grow as supersonic deflagrations. 

Such deflagrations consist of a Jouguet deflagration, followed by a rarefaction 

wave. Depending on the amount of supercooling, the maximal velocity of 

supersonic deflagrations varies between the sound and the light velocities. 

The solutions faster than supersonic deflagrations are weak detonations. 
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I. INTRODUCTION 



When an exothermic process takes place in an ideal fluid, knowing the hydrodynamics of 
growing bubbles becomes relevant. Such situations include spherically symmetric chemical 
burning in combustible gases , and the growth of bubbles of the low-temperature phase 
in cosmological flrst-order phase transitions P-P^. The solutions of the energy-momentum 
conservation equations across the 'burning' surface can naturally be divided into six hydro- 
dynamical modes with different qualitative features, called the weak, strong, and Jouguet 
deflagrations, and the weak, strong, and Jouguet detonations. When stabihty considera- 
tions, the microscopic mechanism of entropy production at the burning surface, and the 
overall boundary conditions, are taken into account, it turns out that some of the hydrody- 
namical modes are actually forbidden. It is the purpose of the present paper to complete 
the existing classiflcation of the allowed hydrodynamical modes in cosmological phase tran- 
sitions, and to present, to our knowledge, a new flow proflle by which the burning surface 
can be connected to the boundary conditions. 

Our new results concern deflagrations, which are the slower of the two main types of 
solutions. Usually, it is assumed that the flow proflle of a deflagration bubble is such that the 
burning surface leaves the matter at rest. Then weak deflagrations are subsonic and Jouguet 
deflagrations expand at sound velocity. The strong deflagrations would be supersonic, but it 
is known for chemical burning that such solutions do not exist . We flnd that the strong 
deflagrations are forbidden in cosmological phase transitions, as well. However, there is 
another type of a supersonic mode. Indeed, a Jouguet deflagration can also leave behind it a 
rarefaction wave, in which case matter is not at rest immediately behind the burning surface. 
Such solutions are supersonic, and flll in the velocity gap between the sound velocity and 
the slowest detonations. As far as we can see, supersonic deflagrations should be possible 
both in chemical burning, and in cosmological phase transitions, such as the QCD and the 
electroweak phase transition. 

The plan of the paper is the following. In Sec. 0, the properties of bubble growth 
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in relativistic ideal fluid are reviewed both for deflagrations and detonations. In Sec. |T|, 
we use a simple dynamical model to study which hydrodynamical growth mechanisms can 
actually be reached from an initial critical bubble configuration. In particular, supersonic 
deflagrations are found to be possible. The structure of spherical supersonic deflagration 
bubbles is analyzed in detail in Sec. |V|, for the bag equation of state. The conclusions are 
in Sec. 0. 



II. DEFLAGRATIONS AND DETONATIONS 

Let us go into the rest frame of the expanding phase transition surface. Imposing energy- 
momentum conservation across the discontinuity gives the equations 

WqllVq = Whllvh (1) 
WqlqVg +Pq = Wh'jlvl + Ph , (2) 

where the subscripts q and h refer to the inflowing high-temperature ('quark') phase and 
the out-flowing low-temperature ('hadron') phase, respectively {w is the proper enthalpy 
density). We choose the signs so that the velocities are positive, and denote by F the 
quantity Wq'jqVq. It then follows that 

F = > . (3) 

This equation tells us that there are two qualitatively different kinds of solutions. For 
deflagrations, the pressure decreases and the velocity increases across the surface: ph < Pq 
and Vh > Vq. For detonations, the opposite is true. For deflagrations, the inflow velocity Vq 
is subsonic, and for detonations, it is supersonic [[l|,p!4|. Note that for deflagrations, the sign 
of the pressure difference is opposite to what it would have been, had the temperature not 
changed across the phase transition surface. 

A further division of the different hydrodynamical processes can be made according to the 
outflow velocity Vh- Both for deflagrations and detonations, energy-momentum conservation 
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allows Vh to be on either side of, or equal to, the sound velocity Cg. If Vh equals c^, the 
process is called a Jouguet process. If Vh is on the same side of Cg as f^, the process is called 
weak] otherwise it is called strong. For instance, for weak deflagrations, Vh < c^. 

To construct a complete bubble solution, the discontinuity at the phase transition surface 
has to be matched with possible shock discontinuities and continuous flow profiles, to satisfy 
the boundary conditions. A deflagration solution can be constructed from a supersonic shock 
front preceding the phase transition surface and setting the cosmic matter moving, followed 
by a continuous flow profile, and terminated by the phase transition front bringing the matter 
back at rest. A detonation solution can be constructed from a weak or Jouguet detonation 
front moving supersonically into matter at rest, and followed by a continuous rarefaction 
wave bringing the matter again at rest. With a strong detonation, the boundary conditions 
of a spherically expanding bubble cannot be satisfied Flow profiles of deflagration 

bubbles are given in Ref. 0, and flow profiles of detonation bubbles are illustrated, e.g., 
in Fig. 2 of Ref. These simplest complete flow profiles are dubbed, according to the 



nature of the phase transition front, the weak, Jouguet and strong deflagrations, and the 
weak and Jouguet detonations. In the next Section, it is seen that the presented solutions 
do not quite exhaust the flow profiles relevant for cosmological phase transitions. In the 
remainder of this Section, on the other hand, it is argued that some of the solutions will not 
be realized in nature. 

First, strong deflagrations are expected to be forbidden. For chemical combustion, there 
are at least two ways of seeing this. Assuming that the burning front can be described as a 
changing mixture of the burnt and the unburnt matter, one can see that strong deflagrations 



would include regions of negative entropy production and they are therefore impos- 

sible. This proof is not valid for cosmological phase transitions, since the assumption of a 
changing mixture would lead to a vanishing surface tension at T^.. There is another way of 
excluding the strong deflagrations, coming from the observation that even if such a solution 
could momentarily be forced to exist, it would be unstable and tend to split into other 
solutions This argument should hold also for cosmological phase transitions. Hence we 



expect to see only weak and Jouguet deflagrations. These expand at most at sound velocity. 

Second, consider weak detonations. To study them, one must know something about 
the microscopic mechanism of phase separation. For chemical combustion the rate 

of burning rises dramatically with temperature, and the velocity of the burning front is 
determined by how fast heat can be supplied to yet unburnt regions. As a consequence, the 
microscopic structure of a detonation front is such that a shock front heats the matter so 
much that rapid burning can begin, and a layer where the burning takes place follows just 
behind the shock frontQ. For weak detonations, the layer behind the shock front is a strong 
deflagration. Hence, the impossibility of strong deflagrations implies the impossibility of 
weak detonations. For cosmological phase transitions, the structure of the burning front is 
different. A phase transition is more apt to start at a low temperature, when supercooling 
is larger. In particular, a point in space outside the original critical bubble does not wait for 
a rise in temperature to start transforming into the new phase, but the transformation is 
caused by interactions with the new phase in the detonation front. Hence, the structure of a 
detonation front need not be a shock front followed by a deflagration, and weak detonations 
cannot be excluded in cosmological phase transitions . 

Let us note that there exist processes in classical hydrodynamics which have the same 
classification of the different hydrodynamical growth mechanisms as cosmological phase 
transitions, unlike chemical burning. Such are condensation discontinuities §132]. An 
example of a condensation discontinuity is the condensation of supersatured water vapour 



^This assumption was also made by Steinhardt in Ref. Q. In Ref. Q weak detonations were ex- 
cluded from consideration by the statement "The lower point b at which the chord ad intersects the 
detonation adiabat cannot be reached if the reaction is exothermic" , which assumes the detonation 
is initiated by a shock. Ref. Q is often quoted as a proof that Jouguet detonations are the only 
kinds of detonations that occur, but this is a result of relativistic combustion theory, which does 
not apply to cosmological phase transitions. 
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into small water droplets, i.e., fog. This is a phase transition, though with a finite chemical 
potential. Indeed, only strong deflagrations and detonations are excluded in Ref. and 
the other processes are accepted. 



III. NUMERICAL RESULTS 

To be able to check the presented arguments about the possibility of the different kinds 
of solutions, one would like to have a way of taking automatically into account energy- 
momentum conservation, boundary conditions, and the microscopic structure of the phase 
transition surface. One might also wish to follow the real-time history of the growing bubbles, 
to see which stationary final states can actually be reached from reasonable initial conditions. 
These problems can be solved by imposing an equation of motion for some continuously 
changing order parameter field (f{t,x.), and by then solving the equation. The equation 
could, in principle, be derived from non-equilibrium field theory. 

In Ref. |]T5[, the simple phenomenological form 

dV 

d^d^if + ^ + W'^d^if = (4) 



of the equation for the yj-field was proposed (in Ref. ITS the notation l/F was used for 



7]). This is the simplest relativistic generalization of the Langevin equation without the 
noise term. The parameter r] fixes the one degree of freedom left free by energy-momentum 
conservation equations. From eq. @ and energy-momentum conservation it follows that the 
entropy density s satisfies 

Td,{su^) = r]{u^d,vf , (5) 

so that 1] can be interpreted as a parametrization of the entropy production at the phase 
transition surface. Below, rj is treated as a free parameter, although for the EW phase transi- 
tion, its value could be estimated by comparing the analytic formula for the defiagration 
front velocity in the limit of small supercooling with microscopic calculations [|T^-|T3 of the 
same quantity. 



In Ref. |]T5[, eq. (^) and the energy- momentum conservation equations were solved nu- 
merically for the planar symmetric case, corresponding to a 1+1 dimensional spacetime. 
Note that the dimension of space d does not affect the results qualitatively, since the phase 
transition front looks locally planar, and the qualitative features of the continuous flow pro- 
files outside the phase transition front do not depend on c? There are two methods of 
solving eq. @). Either one can take some initial bubble, and integrate all the partial differ- 
ential equations forward in time (the dynamical code) , or one can search only for stationary 
final states (the static code). 

The general structure of the stationary final states is shown in Fig. In this figure, 
three kinds of solutions are displayed. First, the dotted lines (a, . . .) each represent a one- 
parameter family of solutions, with a given initial supercooling Tj, which satisfy the energy- 
momentum conservation equations and the boundary conditions, alone. In drawing these 
curves, it was assumed that the solutions are of the type discussed in Sec. 0, Second, the solid 
lines represent the solutions which satisfy energy-momentum conservation and eq. (^ across 
the phase transition surface, with a given 77, but for which the global boundary conditions 
have not been imposed. For given Ty and r], the possible stationary final states are then 
found from the intersections of these two kinds of curves. Third, the gray lines indicate those 
stationary final states which are actually realized by the dynamical code, starting from a 
reasonable initial state, for given Tj, as a function of rj. As illustrated with the arrows for 
the curve a, the solution is a weak deflagration in the limit 77 — > 00, moves towards strong 
deflagrations as 77 is decreased, and changes to a detonation at some point. 

A few interesting things can be seen in Fig. |l[ First, one should note the general feature 
that weak detonations exist, as anticipated in Sec. 0. Second, the dotted lines have a very 
large gradient for deflagrations. This indicates that the temperature Tq in front of the phase 
transition surface may be much higher than Tj due to the heating caused by the shock front, 
but the temperature Th behind the phase transition surface is very near Tf. For detonations, 
similarly, it turns out that the temperature in the center of the expanding bubble differs 
much less from Tf than the temperature just behind the detonation front. Third, it 



is seen that for some Tj and rj, there are more than one possible stationary final states. 
For instance, for rj = O.OlOTc and Ty = 0.98Tc, there are three possible final states. The 
dynamical code selects one of theseQ. The solutions which are not realized have a sharp peak 
in temperature expanding at roughly the speed of sound. If one gives the unreached weak 
defiagration as an initial condition to the dynamical code, it remains as such, indicating 
that this solution may be reachable from other initial conditions. However, the unreached 
weak detonation transforms to either of the other solutions, suggesting an instability. 

Let us now turn to the existence of supersonic defiagrations. From Fig. |l], it is seen that 
strong defiagrations have not been found with the static code. Technically, this is due to the 
fact that no matter how well one tries to guess the correct solution, the guess does not relax 
to a strong defiagration but rather changes considerably to form a weak detonation, with the 
same T^. This is in accordance with strong defiagrations being unstable, as argued in the 
previous Section. On the other hand, with the dynamical code, we have found supersonic 
solutions which lie very near the dotted lines representing strong defiagrations (see curve 
a). A careful inspection reveals, however, that in the region covered by the lighter shade, 
the gray curve is not exactly on the dotted curve. By drawing a figure of the temperature 
profiles, the reason becomes apparent. In Fig. ^ temperature profiles are shown for the 
supercooling Tf = O.SOTc. The supersonic defiagration solution {rj = 0.17Tc) is seen to be 
comprised of a shock front, a phase transition front, and a rarefaction wave in the low- 
temperature phase. If the temperature just behind the phase transition front is taken as the 
x-coordinate in Fig. |l] instead of the temperature in the center of the bubble, the solution 
lies on the dashed curve corresponding to a Jouguet defiagration. Furthermore, if one gives 
a strong defiagration as an initial condition to the dynamical code, the solution is unstable, 
and transforms into a Jouguet defiagration followed by a rarefaction wave. In the next 



^Figs. 19 and 20 of Ref. represent one solution of the static code. The solution of the dynamical 
code actually changes into a weak detonation at around rj = 0.03Tc, with Vdct ~ 0.7. 
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Section, these solutions are analyzed in more detail in 1+3 dimensions. 



IV. SPHERICAL SOLUTIONS 



For an overview of the different hydrodynamical solutions for spherical phase transition 
bubbles, we employ a simpler equation of state, the bag EOS 



where B is the bag constant (e is the proper energy density). We assume a spherically 
symmetric similarity solution, i.e., it can be described in terms of the single coordinate 
^ = r/t. In contrast to the previous Section, we are here not interested in how the bubble 
relaxes to the similarity solution from some initial condition. We also give up modelling the 
microscopic structure of the phase transition front, and treat it as a discontinuity. Therefore 
we do not have a parameter rj coming from microphysics to select among the different 
solutions hydrodynamically allowed. We just parametrize the solutions by the velocity of 
the phase transition front ^pt, (,^dcfi or ^dct)- Then we get almost analytical {i.e., involving 
ordinary differential equations) solutions for the bubble structure, instead of having to resort 
to numerical hydrodynamics. 

The continuous parts of the flow satisfy the equations p 



1 

Ph = -eh, 



(6) 




(7) 



dv 



2 



(8) 



d^ 



and 



1 de 
e d^ 



4 



1-^vl-v^dC 



^ — V 1 dv 



(9) 



where e = e in the h phase and e = e — B in the q phase. 



In Fig. ^ we show a sequence of profiles for phase transition bubbles, with the same 
initial energy density ej, but for different ^p t.. 
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If the phase transition front is subsonic, ^p.t. = ^defi < c^, we have a deflagration bubble 
of the type described in 0. The phase transition front is preceded by a shock at > c^, 
the matter is at rest at ,^ > ^sh and at ,^ < ^defl, and the shocked fluid in between has a 
continuous flow profile, with the (outward) flow velocity and energy density increasing with 
decreasing ^. We denote by ^sh,j the shock velocity of the Jouguet deflagration bubble with 
^p.t. = Cg. The slowest deflagrations have extremely weak shocks at ^sh ~ Cg. 

The slowest detonation is the Jouguet detonation, with 



Solutions with t. > ^det,j are weak detonations. 

The velocity range between ^p.t. = Cs and ^p.t. = ^det,j is covered by the supersonic 
deflagrations. As ^p.t. is increased past Cs, a rarefaction wave appears in the region between 
i = Cs and ^ = ^p.t.. The phase transition remains a Jouguet deflagration, and is preceded by 
a shock, with ^p t. < .^sh < ■Cdct.j- As ,^p.t. is increased, ^sh grows also, but not as rapidly, and 
thus the region of shocked fluid becomes narrower. The two velocities ^p,t. and become 
equal at ^p.t. = ^sh = ^det,j, the region of shocked fluid disappears and the solution becomes 
a Jouguet detonation. 

In the limit where ^p.t. approaches ^det,j from below, the shocked region becomes infinites- 
imally thin. This resembles the microscopic structure of a detonation front in chemical com- 
bustion, where heating by a shock initiates the transition. However, the weak detonations 
(^p.t. > ^det,j) here do not have such structure. 

In Fig. ^ we show the structure of a supersonic deflagration bubble. Figs. ^ and |^ show 
how the different quantities in the bubble change as a function of ^p.t.. Fig. ^ shows the 
limiting velocities ^sh,j and ^det,j as a function of ^f/B. 

Let us note that in this Section we have studied pure hydrodynamics, whereas one 
should actually also insist on entropy production being non-negative at the phase transition 
front 1^. This requires knowledge of the dimensionless parameter r = Og/a/j in the bag EOS 
(e/i = Sa/iT"^, eg = ?>aqT^+B). The parameter r fixes the energy density e/i? = (4r— l)/(r— 1) 




(10) 
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corresponding to the critical temperature Tc- The entropy condition restricts the availabihty 
of the different hydrodynamical modes near T^, and for a sufficiently small supercooling, only 
weak deflagrations remain allowed ||I9| , Sec. IV]. However, considering any e/, for a sufficiently 
small r this value of e/ will correspond to a large supercooling, and all the hydrodynamical 
modes are allowed. Therefore the entropy condition does not change our conclusions. 

V. CONCLUSIONS 

We have discovered a new kind of a bubble solution for cosmological phase transitions. 
The phase transition front of such a bubble is a Jouguet deflagration, is preceded by a shock, 
and followed by a rarefaction wave. The phase transition front moves supersonically with 
respect to the matter at rest. These supersonic deflagration bubbles fill the velocity gap 
between weak deflagrations, which are subsonic, and weak detonations. The two limiting 
cases of the supersonic deflagration bubble are the 'ordinary' Jouguet deflagration bubble 
[i.e., without a rarefaction wave) and the Jouguet detonation bubble. 

It has been known that strong deflagrations would not occur since they are unstable to 
1-dimensional perturbations {i.e., in the radial direction). By following a strong deflagra- 
tion bubble with a numerical hydrodynamics code we have seen how such a perturbation 
transforms the strong deflagration into a Jouguet deflagration by developing a rarefaction 
wave. This gives rise to the supersonic deflagration bubble we have described. 

By modelling the microscopic mechanism of entropy production at the phase transition 
front by a simple dynamical model, we found that the solutions 'near' Jouguet detonations 
[e.g., proflles 7-9 in Fig. Q were in many cases not realized in a numerical evolution, al- 
though they were allowed by the external conditions (including the entropy condition). The 
velocity ^p.t. of the phase transition front was found not to increase continuously as the 
microscopic entropy-production parameter rj was decreased, but the solution jumped from a 
weak deflagration to a much faster weak detonation. Thus the most 'violent' solutions, i.e., 
the ones with large energy densities and flow velocities occurring near the phase transition 
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front, were skipped over. This happened in the cases with a small supercooling, and we had 
to go to large supercooling to see examples of supersonic deflagrations. Examining Fig. |1|, 
one finds for small supercooling that, for a given value of Tf and rj there can be as many as 
three solutions allowed, not just one. One of these is a deflagration, and two are detonations. 
In the numerical evolution, the weaker, i.e., the faster, of the detonations was then realized. 

This suggests that in nature the more dramatic solutions, the Jouguet detonation and 
solutions near it, are more difficult to realize, especially for small supercooling. This is in 
sharp contrast to chemical burning, where the Jouguet detonations are the only possible 
detonation solutions. 
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FIGURES 



FIG. 1. The solutions of the model of Ref. |15], in the case of planar symmetry. In 1+3 
dimensions, the solid lines remain the same, but the exact functional form of the dotted lines 
changes a little. For detonations, the x-axis is the temperature just behind the phase transition 
surface; for deflagrations, it is the temperature in the center of the expanding bubble. The y-axis is 
the temperature just in front of the phase transition surface. The region forbidden kinematically, 
by the non-negativity of entropy production, or by boundary conditions, has been covered with 
the darker shade, assuming the types of solutions discussed in Sec. |IJ The two allowed regions 
represent the deflagrations and the weak detonations. Dashed lines indicate the Jouguet processes, 
and strong deflagrations are covered with the lighter shade. The meaning of the different curves 
has been explained in the text. The nucleation temperatures Tf/Tc corresponding to a, /?, 7, 
5 and e are 0.86, 0.89, 0.92, 0.95 and 0.98. The (QCD-type) parameters used in this figure are 
L = O.IT^, a = O.IT^ and Ic = ST"^. Given the expansion rate of the Universe, the actual 
nucleation temperature Tj/Tc is then about 0.86 (Fig. 1 in Ref. ||l9|). 

FIG. 2. Temperature profiles at t = 14400r~^ for different values of rj. All the profiles are 
moving to the right. The arrows indicate the location of the phase transition front. For rj = 0.20Tc, 
the solution is a weak deflagration; for rj = O.lTTc, it is a Jouguet deflagration followed by a 
rarefaction wave; for -q = 0.16rc and rj = 0.12Tc, the solution is a weak detonation. The velocity 
of the phase transition front changes from 0.56 at 77 = 0.20Tc to 0.74 at r/ = 0.12Tc. 

FIG. 3. Energy-density (upper set) and flow velocity (lower set) profiles for spherical phase 
transition bubbles. We show 13 profiles, all for the same initial energy density £/ = lOOi?, but with 
different phase transition front velocities. The six slowest (on the left) are weak deflagrations. The 
seventh is a Jouguet deflagration bubble, with ^dcfi = Cg. The eighth is a supersonic deflagration 
bubble, which has a Jouguet deflagration front preceded by a shock and followed by a rarefaction 
wave. The ninth is a supersonic deflagration bubble almost at the limit where it becomes a Jouguet 
detonation. The shocked region is now so thin that it appears as a mere vertical line in the figure. 
The last four (on the right) are weak detonations. 
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FIG. 4. A supersonic deflagration bubble. The upper curve is the energy density profile and the 
lower curve the flow velocity proflle. This flgure is for a supersonic spherical deflagration bubble 
with an initial energy density £/ = lOOB and shock velocity ^gh = 0.625. The velocity of the phase 
transition front is ^defl = 0.6163. 

FIG. 5. The energy densities (a) and flow velocities (b) at special points of the bubble profile 
(see Fig. 4). These quantities are shown as a function of the phase transition velocity ^p.t. for a 
fixed initial energy density £/ = 100-B. The solution switches from a weak deflagration bubble 
to a supersonic deflagration bubble at ^p.t. = Cg = 0.57735 and to a weak detonation bubble at 
^p.t. = Cdet.J = 0.6534. 

FIG. 6. The special points of the bubble profile as a function of ^p.t. for e/ = lOOB. A weak 
deflagration bubble has a region of shocked fluid between ^ = ^defl and ^ = ^gh- A weak detonation 
bubble has a rarefaction wave extending from ^ = to ^ = ^det • A supersonic deflagration bubble 
has both features. 

FIG. 7. The velocity ^det,J of the phase transition front in the Jouguet detonation bubble 
and the shock velocity ^sh,j in the Jouguet deflagration bubble as a function of the initial energy 
density e/. As the velocity of a weak deflagration is increased from to c^, the shock velocity 
increases from Cg to ^sh,j, at which point the solution becomes a supersonic deflagration bubble. 
As the deflagration velocity is further increased from Cg to ^det,J the shock velocity increases from 
Csh,j to ^det,J at which point the phase transition front catches up with the shock and the solution 
becomes a detonation bubble. 
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The classification of the hydrodynamical growth mecha- 
nisms for the spherical bubbles of the low-temperature phase 
in cosmological phase transitions is completed by showing 
that the bubbles can grow as supersonic deflagrations. Such 
deflagrations consist of a Jouguet deflagration, followed by a 
rarefaction wave. Depending on the amount of supercooling, 
the maximal velocity of supersonic deflagrations varies be- 
tween the sound and the light velocities. The solutions faster 
than supersonic deflagrations are weak detonations. 
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I. INTRODUCTION 

When an exothermic process takes place in an ideal 
fluid, knowing the hydrodynamics of growing bubbles be- 
comes relevant. Such situations include spherically sym- 
metric chemical burning in combustible gases [1,2], and 
the growth of bubbles of the low-temperature phase in 
cosmological first-order phase transitions [3-13]. The so- 
lutions of the energy-momentum conservation equations 
across the 'burning' surface can naturally be divided into 
six hydrodynamical modes with different qualitative fea- 
tures, called the weak, strong, and Jouguet deflagrations, 
and the weak, strong, and Jouguet detonations. When 
stability considerations, the microscopic mechanism of 
entropy production at the burning surface, and the over- 
all boundary conditions, are taken into account, it turns 
out that some of the hydrodynamical modes are actu- 
ally forbidden. It is the purpose of the present paper to 
complete the existing classification of the allowed hydro- 
dynamical modes in cosmological phase transitions, and 
to present, to our knowledge, a new flow profile by which 
the burning surface can be connected to the boundary 
conditions. 

Our new results concern deflagrations, which are the 
slower of the two main types of solutions. Usually, it is 
assumed that the flow profile of a deflagration bubble is 
such that the burning surface leaves the matter at rest. 
Then weak deflagrations are subsonic and Jouguet defla- 
grations expand at sound velocity. The strong deflagra- 
tions would be supersonic, but it is known for chemical 
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burning that such solutions do not exist [1,2]. We find 
that the strong deflagrations are forbidden in cosmologi- 
cal phase transitions, as well. However, there is another 
type of a supersonic mode. Indeed, a Jouguet deflagra- 
tion can also leave behind it a rarefaction wave, in which 
case matter is not at rest immediately behind the burn- 
ing surface. Such solutions are supersonic, and fill in the 
velocity gap between the sound velocity and the slowest 
detonations. As far as we can see, supersonic deflagra- 
tions should be possible both in chemical burning, and 
in cosmological phase transitions, such as the QCD and 
the electroweak phase transition. 

The plan of the paper is the following. In Sec. II, 
the properties of bubble growth in relativistic ideal fluid 
are reviewed both for deflagrations and detonations. In 
Sec. Ill, we use a simple dynamical model to study 
which hydrodynamical growth mechanisms can actually 
be reached from an initial critical bubble configuration. 
In particular, supersonic deflagrations are found to be 
possible. The structure of spherical supersonic deflagra- 
tion bubbles is analyzed in detail in Sec. IV, for the bag 
equation of state. The conclusions are in Sec. V. 

II. DEFLAGRATIONS AND DETONATIONS 

Let us go into the rest frame of the expanding phase 
transition surface. Imposing energy-momentum conser- 
vation across the discontinuity gives the equations 

'Wq'Tq'Vq = Whflvh (1) 
■^975^5 +Vq= WhlWh + Vh , (2) 

where the subscripts q and h refer to the inflowing 
high-temperature ('quark') phase and the out-flowing 
low-temperature ('hadron') phase, respectively (w is the 
proper enthalpy density). We choose the signs so that 
the velocities are positive, and denote by F the quantity 
Wqj^Vq. It then follows that 

F = > . (3) 

Vq - Vh 

This equation tells us that there are two qualitatively dif- 
ferent kinds of solutions. For deflagrations, the pressure 
decreases and the velocity increases across the surface: 
Ph < Pq and Vfi > Vq. For detonations, the opposite is 
true. For deflagrations, the inflow velocity Vq is subsonic, 
and for detonations, it is supersonic [1,14]. Note that for 
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deflagrations, the sign of the pressure difference is oppo- 
site to what it would have been, had the temperature not 
changed across the phase transition surface. 

A further division of the different hydrodynamical 
processes can be made according to the outflow veloc- 
ity Vh. Both for deflagrations and detonations, energy- 
momentum conservation allows Vh to be on either side 
of, or equal to, the sound velocity Cj. livh equals Cj, the 
process is called a Jouguet process. If Vh is on the same 
side of Cs as Vq, the process is called weak; otherwise 
it is called strong. For instance, for weak deflagrations. 

To construct a complete bubble solution, the disconti- 
nuity at the phase transition surface has to be matched 
with possible shock discontinuities and continuous flow 
profiles, to satisfy the boundary conditions. A defla- 
gration solution can be constructed from a supersonic 
shock front preceding the phase transition surface and 
setting the cosmic matter moving, followed by a continu- 
ous flow profile, and terminated by the phase transition 
front bringing the matter back at rest. A detonation 
solution can be constructed from a weak or Jouguet det- 
onation front moving supersonically into matter at rest, 
and followed by a continuous rarefaction wave bringing 
the matter again at rest. With a strong detonation, the 
boundary conditions of a spherically expanding bubble 
cannot be satisfied [1,3]. Flow profiles of deflagration 
bubbles are given in Ref. [5] , and flow profiles of detona- 
tion bubbles are illustrated, e.g., in Fig. 2 of Ref. [14]. 
These simplest complete flow profiles are dubbed, ac- 
cording to the nature of the phase transition front, the 
weak, Jouguet and strong deflagrations, and the weak 
and Jouguet detonations. In the next Section, it is seen 
that the presented solutions do not quite exhaust the 
flow profiles relevant for cosmological phase transitions. 
In the remainder of this Section, on the other hand, it is 
argued that some of the solutions will not be realized in 
nature. 

First, strong deflagrations are expected to be forbid- 
den. For chemical combustion, there are at least two 
ways of seeing this. Assuming that the burning front 
can be described as a changing mixture of the burnt and 
the unburnt matter, one can see that strong deflagra- 
tions would include regions of negative entropy produc- 
tion [2,14], and they are therefore impossible. This proof 
is not valid for cosmological phase transitions, since the 
assumption of a changing mixture would lead to a van- 
ishing surface tension at Tc. There is another way of 
excluding the strong deflagrations, coming from the ob- 
servation that even if such a solution could momentarily 
be forced to exist, it would be unstable and tend to split 
into other solutions [1]. This argument should hold also 
for cosmological phase transitions. Hence we expect to 
see only weak and Jouguet deflagrations. These expand 
at most at sound velocity. 

Second, consider weak detonations. To study them, 
one must know something about the microscopic mecha- 
nism of phase separation. For chemical combustion [1,2], 



the rate of burning rises dramatically with temperature, 
and the velocity of the burning front is determined by 
how fast heat can be supplied to yet unburnt regions. 
As a consequence, the microscopic structure of a detona- 
tion front is such that a shock front heats the matter so 
much that rapid burning can begin, and a layer where the 
burning takes place follows just behind the shock front^. 
For weak detonations, the layer behind the shock front is 
a strong deflagration. Hence, the impossibility of strong 
deflagrations implies the impossibility of weak detona- 
tions. For cosmological phase transitions, the structure 
of the burning front is different. A phase transition is 
more apt to start at a low temperature, when supercool- 
ing is larger. In particular, a point in space outside the 
original critical bubble does not wait for a rise in tem- 
perature to start transforming into the new phase, but 
the transformation is caused by interactions with the new 
phase in the detonation front. Hence, the structure of a 
detonation front need not be a shock front followed by a 
deflagration, and weak detonations cannot be excluded 
in cosmological phase transitions [14]. 

Let us note that there exist processes in classical hy- 
drodynamics which have the same classification of the 
different hydrodynamical growth mechanisms as cosmo- 
logical phase transitions, unlike chemical burning. Such 
are condensation discontinuities [1, §132]. An example of 
a condensation discontinuity is the condensation of super- 
satured water vapour into small water droplets, i.e., fog. 
This is a phase transition, though with a finite chemical 
potential. Indeed, only strong deflagrations and detona- 
tions are excluded in Ref. [1], and the other processes are 
accepted. 

III. NUMERICAL RESULTS 

To be able to check the presented arguments about 
the possibility of the different kinds of solutions, one 
would like to have a way of taking automatically into 
account energy-momentum conservation, boundary con- 
ditions, and the microscopic structure of the phase transi- 
tion surface. One might also wish to follow the real-time 
history of the growing bubbles, to see which stationary fi- 
nal states can actually be reached from reasonable initial 
conditions. These problems can be solved by imposing 
an equation of motion for some continuously changing 
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In Ref. [3] weak detonations were excluded from considera- 
tion by the statement "The lower point b at which the chord 
ad intersects the detonation adiabat cannot be reached if the 
reaction is exothermic", which assumes the detonation is ini- 
tiated by a shock. Ref. [3] is often quoted as a proof that 
Jouguet detonations are the only kinds of detonations that 
occur, but this is a result of relativistic combustion theory, 
which does not apply to cosmological phase transitions. 
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order parameter field ^(i,x), and by then solving the 
equation. The equation could, in principle, be derived 
from non-equilibrium field theory. 

In Ref. [15], the simple phenomenological form 

dV 

d^d^ip + -^ + Vn^d^v = (4) 

of the equation for the ^-field was proposed (in Ref. [15] 
the notation l/F was used for 77). This is the simplest 
relativistic generalization of the Langevin equation with- 
out the noise term. The parameter 77 fixes the one degree 
of freedom left free by energy-momentum conservation 
equations. From eq. (4) and energy-momentum conser- 
vation it follows that the entropy density s satisfies 

Td^isu'') = v{n^d^vf , (5) 

so that 77 can be interpreted as a parametrization of the 
entropy production at the phase transition surface. Be- 
low, 77 is treated as a free parameter, although for the 
EW phase transition, its value could be estimated [15] 
by comparing the analytic formula for the deflagration 
front velocity in the limit of small supercooling with mi- 
croscopic calculations [16-18] of the same quantity. 

In Ref. [15], eq. (4) and the energy-momentum con- 
servation equations were solved numerically for the pla- 
nar symmetric case, corresponding to a 1+1 dimensional 
spacetime. Note that the dimension of space d does not 
affect the results qualitatively, since the phase transition 
front looks locally planar, and the qualitative features of 
the continuous flow profiles outside the phase transition 
front do not depend on d [5]. There are two methods of 
solving eq. (4). Either one can take some initial bubble, 
and integrate all the partial differential equations forward 
in time (the dynamical code), or one can search only for 
stationary final states (the static code). 

The general structure of the stationary final states is 
shown in Fig. 1. In this figure, three kinds of solutions are 
displayed. First, the dotted lines (a, . . .) each represent a 
one-parameter family of solutions, with a given initial su- 
percooling Tf , which satisfy the energy-momentum con- 
servation equations and the boundary conditions, alone. 
In drawing these curves, it was assumed that the solu- 
tions are of the type discussed in Sec. II, Second, the 
solid lines represent the solutions which satisfy energy- 
momentum conservation and eq. (4) across the phase 
transition surface, with a given 77, but for which the 
global boundary conditions have not been imposed. For 
given Tf and 77, the possible stationary final states are 
then found from the intersections of these two kinds of 
curves. Third, the gray lines indicate those stationary 
final states which are actually realized by the dynamical 
code, starting from a reasonable initial state, for given 
Ty , as a function of 77. As illustrated with the arrows for 
the curve a, the solution is a weak deflagration in the 
limit 77 ^ 00, moves towards strong deflagrations as 77 is 
decreased, and changes to a detonation at some point. 




FIG. 1. The solutions of the model of Ref. [15], in the case 
of planar symmetry. In 1+3 dimensions, the solid lines re- 
main the same, but the exact functional form of the dotted 
lines changes a little. For detonations, the E-axis is the tem- 
perature just behind the phase transition surface; for defla- 
grations, it is the temperature in the center of the expanding 
bubble. The 2/-axis is the temperature just in front of the 
phase transition surface. The region forbidden kinematically, 
by the non-negativity of entropy production, or by boundary 
conditions, has been covered with the darker shade, assuming 
the types of solutions discussed in Sec. 11. The two allowed 
regions represent the deflagrations and the weak detonations. 
Dashed lines indicate the Jouguet processes, and strong de- 
flagrations are covered with the lighter shade. The meaning 
of the different curves has been explained in the text. The 
nucleation temperatures TfjTc corresponding to a, /3, 7, 6 
and e are 0.86, 0.89, 0.92, 0.95 and 0.98. The (QCD-type) 
parameters used in this figure are L = O.IT^, a = O.IT^ and 
Ic = 6Te~^ . Given the expansion rate of the Universe, the ac- 
tual nucleation temperature Tf/Tc is then about 0.86 (Fig. 1 
in Ref. [19]). 

A few interesting things can be seen in Fig. 1. First, 
one should note the general feature that weak detona- 
tions exist, as anticipated in Sec. II. Second, the dotted 
lines have a very large gradient for deflagrations. This 
indicates that the temperature Tq in front of the phase 
transition surface may be much higher than Tf due to 
the heating caused by the shock front, but the tempera- 
ture Tfi behind the phase transition surface is very near 
Tf . For detonations, similarly, it turns out that the tem- 
perature in the center of the expanding bubble differs 
much less from Tf than the temperature Th just behind 
the detonation front. Third, it is seen that for some Tf 
and 77, there are more than one possible stationary final 
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states. For instance, for 77 = O.OlOTc and Tf = 0.98Tc, 
there are three possible final states. The dynamical code 
selects one of these^. The solutions which are not realized 
have a sharp peak in temperature expanding at roughly 
the speed of sound. If one gives the unreached weak de- 
flagration as an initial condition to the dynamical code, 
it remains as such, indicating that this solution may be 
reachable from other initial conditions. However, the un- 
reached weak detonation transforms to either of the other 
solutions, suggesting an instability. 

Let us now turn to the existence of supersonic defla- 
grations. From Fig. 1, it is seen that strong deflagrations 
have not been found with the static code. Technically, 
this is due to the fact that no matter how well one tries 
to guess the correct solution, the guess does not relax to 
a strong deflagration but rather changes considerably to 
form a weak detonation, with the same T^. This is in 
accordance with strong deflagrations being unstable, as 
argued in the previous Section. On the other hand, with 
the dynamical code, we have found supersonic solutions 
which lie very near the dotted lines representing strong 
deflagrations (see curve a). A careful inspection reveals, 
however, that in the region covered by the lighter shade, 
the gray curve is not exactly on the dotted curve. By 
drawing a figure of the temperature profiles, the reason 
becomes apparent. In Fig. 2, temperature profiles are 
shown for the supercooling Tf = O.SOTc. The supersonic 
deflagration solution (77 = O.lTTc) is seen to be comprised 
of a shock front, a phase transition front, and a rarefac- 
tion wave in the low-temperature phase. If the temper- 
ature just behind the phase transition front is taken as 
the s-coordinate in Fig. 1 instead of the temperature in 
the center of the bubble, the solution lies on the dashed 
curve corresponding to a Jouguet deflagration. Further- 
more, if one gives a strong deflagration as an initial con- 
dition to the dynamical code, the solution is unstable, 
and transforms into a Jouguet deflagration followed by a 
rarefaction wave. In the next Section, these solutions are 
analyzed in more detail in 1+3 dimensions. 



^Figs. 19 and 20 of Ref. [15] represent one solution of the 
static code. The solution of the dynamical code actually 



changes into a weak detonation at around 77 
Vi^i w 0.7. 



0.03Te, with 



T / 
0.86 - 



0.84 



0.82 



0.80 



7] = 0.20T, 

ri = 0.17T, 

r] = 0.16T, 

— - T] = 0.12T, 



T- 



8000 



9000 



10000 



X / 



FIG. 2. Temperature profiles at t = 14400Te"^ for different 
values of 77. All the profiles are moving to the right. The 
arrows indicate the location of the phase transition front. For 
77 = 0.20Tc, the solution is a weak deflagration; for 77 = O.lTTc, 
it is a Jouguet deflagration followed by a rarefaction wave; for 
77 = 0.16Tc and 77 = 0.12Tc, the solution is a weak detonation. 
The velocity of the phase transition front changes from 0.56 
at 77 = 0.20Te to 0.74 at 77 = 0.12Te. 



IV. SPHERICAL SOLUTIONS 

For an overview of the different hydrodynamical solu- 
tions for spherical phase transition bubbles, we employ a 
simpler equation of state, the bag EOS 



Ph = -eh, 



^(6,-45), 



(6) 
(7) 



where B is the bag constant (e is the proper energy den- 
sity). We assume a spherically symmetric similarity so- 
lution, i.e., it can be described in terms of the single co- 
ordinate ^ = r/t. In contrast to the previous Section, we 
are here not interested in how the bubble relaxes to the 
similarity solution from some initial condition. We also 
give up modelling the microscopic structure of the phase 
transition front, and treat it as a discontinuity. Therefore 
we do not have a parameter 77 coming from microphysics 
to select among the different solutions hydrodynamically 
allowed. We just parametrize the solutions by the ve- 
locity of the phase transition front ^p.t. (^defl or ^det)- 
Then we get almost analytical (i.e., involving ordinary 
differential equations) solutions for the bubble structure, 
instead of having to resort to numerical hydrodynamics. 

The continuous parts of the flow satisfy the equa- 
tions [5] 
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and 



1 de ^ — V 1 dv 
Idi ~ '^l-^v 1-v^di' 



(8) 



(9) 



where e = e in the h phase and e = e — B in the q phase. 

In Fig. 3 we show a sequence of profiles for phase tran- 
sition bubbles, with the same initial energy density ey, 
but for different ^p.t. . 
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FIG. 3. Energy-density (upper set) and flow velocity (lower 
set) profiles for spherical phase transition bubbles. We show 
13 profiles, all for the same initial energy density e/ = 100-B, 
but with different phase transition front velocities. The six 
slowest (on the left) are weak deflagrations. The seventh is a 
Jouguet deflagration bubble, with ^defl = c^. The eighth is a 
supersonic deflagration bubble, which has a Jouguet deflagra- 
tion front preceded by a shock and followed by a rarefaction 
wave. The ninth is a supersonic deflagration bubble almost at 
the limit where it becomes a Jouguet detonation. The shocked 
region is now so thin that it appears as a mere vertical line in 
the figure. The last four (on the right) are weak detonations. 



^ < ^defl) and the shocked fluid in between has a contin- 
uous flow profile, with the (outward) flow velocity and 
energy density increasing with decreasing ^. We denote 
by ^sh,J the shock velocity of the Jouguet deflagration 
bubble with ^p.t. = Cj. The slowest deflagrations have 
extremely weak shocks at ~ Cs- 

The slowest detonation is the Jouguet detonation, with 
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Solutions with ^p.t. > ^det,J are weak detonations. 

The velocity range between ^p.t. = Cj and ^p.t. = ^det,J 
is covered by the supersonic deflagrations. As ^p.t. is in- 
creased past Cj, a rarefaction wave appears in the region 
between ^ = Cj and ^ = ^p.t.. The phase transition re- 
mains a Jouguet deflagration, and is preceded by a shock, 
with ^p.t. < ^sh < ^det.J- As ^p.t. is increased, ^sh grows 
also, but not as rapidly, and thus the region of shocked 
fluid becomes narrower. The two velocities ^p.t. and 
become equal at ^p.t. = ^sh = ^det.Jj the region of shocked 
fluid disappears and the solution becomes a Jouguet det- 
onation. 

In the limit where ^p.t. approaches ^det,J from below, 
the shocked region becomes infinitesimally thin. This re- 
sembles the microscopic structure of a detonation front 
in chemical combustion, where heating by a shock ini- 
tiates the transition. However, the weak detonations 
(^p.t. > ^det,j) here do not have such structure. 

In Fig. 4 we show the structure of a supersonic de- 
flagration bubble. Figs. 5 and 6 show how the different 
quantities in the bubble change as a function of ^p.t. . 
Fig. 7 shows the limiting velocities ^sh,J and ^det,J as a 
function ofey/B. 

Let us note that in this Section we have studied pure 
hydrodynamics, whereas one should actually also in- 
sist on entropy production being non-negative at the 
phase transition front [4]. This requires knowledge of 
the dimensionless parameter r = aq/ah in the bag EOS 
(eji = 3o/iT'*, eg = 3aqT'^ + B). The parameter r fixes 
the energy density e/B = (4r — l)/(r — 1) corresponding 
to the critical temperature Tc- The entropy condition 
restricts the availability of the different hydrodynamical 
modes near Tc, and for a sufficiently small supercool- 
ing, only weak deflagrations remain allowed [19, Sec. IV]. 
However, considering any ey, for a sufficiently small r this 
value of ey will correspond to a large supercooling, and 
all the hydrodynamical modes are allowed. Therefore the 
entropy condition does not change our conclusions. 



If the phase transition front is subsonic, ^p.t. = ^defl < 
Cs , we have a deflagration bubble of the type described 
in [5]. The phase transition front is preceded by a shock 
at ^sh > Cj, the matter is at rest at ^ > and at 
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FIG. 4. A supersonic deflagration bubble. The upper curve 
is the energy density profile and the lower curve the flow ve- 
locity profile. This figure is for a supersonic spherical defla- 
gration bubble with an initial energy density e/ = lOOB and 
shock velocity ^sh = 0.625. The velocity of the phase transi- 
tion front is ^defl = 0.6163. 
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FIG. 5. The energy densities (a) and flow velocities (b) 
at special points of the bubble profile (see Fig. 4). These 
quantities are shown as a function of the phase transition 
velocity ^p.t. for a fixed initial energy density = lOOB. 
The solution switches from a weak deflagration bubble to a 
supersonic deflagration bubble at ^p.t. = Cs = 0.57735 and to 
a weak detonation bubble at ^p.t. = ^det.j = 0.6534. 
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FIG. 6. The special points of the bubble profile as a func- 
tion of ^p.t. for 6/ = lOOB. A weak deflagration bubble has 
a region of shocked fluid between ^ = ^defl and ^ = ^sh- 
A weak detonation bubble has a rarefaction wave extending 
from ^ = Cs to ^ = ^det- A supersonic deflagration bubble has 
both features. 
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FIG. 7. The velocity ^det.j of the phase transition front in 
the Jouguet detonation bubble and the shock velocity ^sh,j 
in the Jouguet deflagration bubble as a function of the initial 
energy density e/. As the velocity of a weak deflagration is 
increased from to Cs, the shock velocity increases from Cs 
to ^sh,j, at which point the solution becomes a supersonic 
deflagration bubble. As the deflagration velocity is further 
increased from Cs to ^det.j the shock velocity increases from 
^sh,j to ^det.j at which point the phase transition front catches 
up with the shock and the solution becomes a detonation 
bubble. 



V. CONCLUSIONS 

We have discovered a new kind of a bubble solution 
for cosmological phase transitions. The phase transition 
front of such a bubble is a Jouguet deflagration, is pre- 
ceded by a shock, and followed by a rarefaction wave. 
The phase transition front moves supersonically with re- 
spect to the matter at rest. These supersonic deflagration 
bubbles fill the velocity gap between weak deflagrations, 
which are subsonic, and weak detonations. The two lim- 
iting cases of the supersonic deflagration bubble are the 
'ordinary' Jouguet deflagration bubble (i.e., without a 
rarefaction wave) and the Jouguet detonation bubble. 

It has been known that strong deflagrations would not 
occur since they are unstable to 1-dimensional pertur- 
bations (i.e., in the radial direction). By following a 
strong deflagration bubble with a numerical hydrody- 
namics code we have seen how such a perturbation trans- 
forms the strong deflagration into a Jouguet deflagration 
by developing a rarefaction wave. This gives rise to the 
supersonic deflagration bubble we have described. 

By modelling the microscopic mechanism of entropy 
production at the phase transition front by a simple 
dynamical model, we found that the solutions 'near' 
Jouguet detonations [e.jr., profiles 7-9 in Fig. 3] were 
in many cases not realized in a numerical evolution, al- 
though they were allowed by the external conditions (in- 
cluding the entropy condition). The velocity ^p.t. of the 
phase transition front was found not to increase contin- 
uously as the microscopic entropy-production parameter 
77 was decreased, but the solution jumped from a weak 
deflagration to a much faster weak detonation. Thus the 
most 'violent' solutions, i.e., the ones with large energy 
densities and flow velocities occurring near the phase 
transition front, were skipped over. This happened in 
the cases with a small supercooling, and we had to go to 
large supercooling to see examples of supersonic deflagra- 
tions. Examining Fig. 1, one finds for small supercooling 
that, for a given value of Ty and 77 there can be as many 
as three solutions allowed, not just one. One of these is a 
deflagration, and two are detonations. In the numerical 
evolution, the weaker, i.e., the faster, of the detonations 
was then realized. 

This suggests that in nature the more dramatic solu- 
tions, the Jouguet detonation and solutions near it, are 
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more difficult to realize, especially for small supercool- 
ing. This is in sharp contrast to chemical burning, where 
the Jouguet detonations are the only possible detonation 
solutions. 
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